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DEFINITIONS  OF  ATTRIBUTES  FOR  LIMB-SCANNING  OR  LIMB-IMAGING 

REMOTE  SENSING  SYSTEMS 


1.  INTRODUCTION 

This  document  contains  descriptions  of  the  attributes,  or  metadata,  to  be  used  to  access  Hafa 
produced  by  limb-scanning  or  limb-imaging  remote  sensing  systems.  While  this  document  is  specific  to 
the  Remote  Atmospheric  and  Ionospheric  Detection  System  (RAIDS)*’^,  the  types  and  descriptions  of 
attributes  and  the  mathematical  formulas  governing  them  are  general.  Thus  the  document  should  be  useful 
for  subsequent  limb-scanning  missions,  some  of  which  will  be  discussed  below.  The  attributes  are  inputs 
to  a  detailed  or  "program"  catalog  for  a  relational  database  management  system  (RDBMS)  designed  by  the 
Naval  Research  Laboratory  (NRL)  Backgrounds  Data  Center  (BDC).  A  program  catalog  consists  of  tables 
of  parameters  describing  each  scan  and  is  the  basis  for  queries  of  the  system  by  the  user.  The  present 
catalog  contains  a  "scans  data  table"  and  a  "mode  data  table,"  the  union  of  which  specifies  the 
characteristics  of  a  particular  platform  scan  or  instrument  wavelength  scan,  corresponding  to  a  given  mode 
of  instrument  activity.  The  attributes  described  herein  are  entries  to  the  scans  and  modes  data  tables. 

The  inunediate  iq>plication  of  the  catalog  is  to  the  storage  and  access  of  data  produced  by  RAIDS. 
RAIDS  consists  of  eight  instruments  which  provide  spectral  intensities  within  the  ultraviolet,  near-infi:ared, 
and  visible  wavelength  ranges  and  covers  a  tangent  altitude  range  of  50-750  km.  This  will  permit  RAIDS 
to  measure  global  distributions  of  the  major  neutral  and  ionized  species  and  important  minor  constituents 
within  the  thermosphere  and  the  ionosphere.  The  mission  will  nominally  last  for  three  years.  The  host 
platform  was  to  be  a  National  Oceanic  and  Atmospheric  Administration  (NOAA)  Television  Infirared 
Observation  Satellite  (I'lKOS-N  series  weather  satellite)^.  Although  the  unfortunate  loss  of  the  NOAA- 13 
satellite  in  August  of  1993  has  precluded  that  possibility,  we  anticipate  that  RAIDS  will  ride  aboard  a 
satellite  with  a  similar  design  and  orbital  parameters.  As  a  result,  diis  document  will  remain  current. 

A  number  of  other  limb-scanning  systems  which  are  similar  to  RAIDS  in  both  capability  and 
deployment  will  begin  fl3dng  in  calendar  year  1996.  Each  will  again  have  a  nominal  lifetime  of  three  years 
and,  at  a  minimum,  will  cover  a  significant  portion  of  the  extreme  ultraviolet  (EUV)  and  far  ultraviolet 
(FUV)  wavelength  ranges  covered  by  RAIDS.  The  Advanced  Research  and  Global  Observation  Satellite 
(ARGOS)  will  fly  in  mid  1996.  The  ARGOS  mission  is  scientific  in  nature,  including  not  only  limb¬ 
scanning  measurements  comparable  to  RAIDS  at  EUV,  FUV,  and  MUV  (middle  ultraviolet)  wavelengths 
but  also  the  High  Resolution  Airglow/Aurora  Spectroscopy  (HIRAAS)  system  with  much  higher  resolution 
than  RAIDS,  plus  limb-imaging  systems,  as  well.  Beginning  as  early  as  1997,  a  sequence  of  United  States 
Air  Force  (USAF)  Defense  Meteorological  Satellite  Program  (DMSP)  missions  will  fly,  lasting  well  past 
the  year  2010.  The  Special  Sensor  Ultraviolet  Limb-Imager  (SSULI)  aboard  each  DMSP  satellite  will 
provide  the  spectral  data  necessary  for  the  near-real-time  assessment  of  the  state  of  the  diermosphere  and 
ionosphere.  The  RAIDS,  ARGOS,  and  SSULI  missions  and  their  anticipated  follow-ons  will  also  provide 
a  consistent  long-term  data  set  with  which  to  determine  signatures  of  global  change. 

The  generation  of  attributes  occurs  when  raw  data  ("Level  0  data")  arrive  from  the  satellite  data 
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processing  activity.  The  Level  0  data  include  satellite  attitude,  ephemeris,  and  health  data  and  scientific 
instrument  data.  In  the  case  of  the  NOAA  TIROS-N  satellite,  for  example,  transmission  of  the  RAIDS 
Level  0  data  to  ground  processing  activities  would  occur  as  part  of  the  global  area  coverage  (GAC)  data 
stream  .  The  ground  activities  woidd  then  decommutate  the  data,  extracting  the  data  relevant  to  RAIDS, 
creating  Level  0  data  files  of  a  specific  format,  and  transmitting  the  files  to  NRL  over  Internet.  We  assume 
that  the  formatting  conventions  used  by  NCAA"  will  hold  for  the  actual  data  handling  activity:  the  level  0 
data  file  consists  of  blocks  of  several  records.  The  first  record  of  a  block  is  the  "General  Housekeeping 
Record"  which  contains  all  satellite  related  data.  The  remaining  records  in  the  block  contain  scientific 
instrument  health  and  observations. 

The  RAIDS  Date  Processing  Facility  (DPF)  converts  the  raw  compressed  Level  0  Hatg  to  the 
archival  Level  1  format  in  several  steps.  The  Level  1  date  set  constitutes  the  primary  archive  and  consists 
of  expanded,  uncalibrated  date.  The  DPF  software  first  expands  the  705  byte  GAC  record  into  844  byte 
records  and  breaks  the  NOAA  files  into  orbit  size  files.  Next  the  software  calculates  the  program  catalog 
attribute  information  and  performs  instrument  health  checks.  The  attribute  information  is  then  given  to  the 
database  system  operator  for  insertion  into  the  Database  Management  System  (DBMS). 


The  following  coordinate  systems  are  used  in  various  attribute  calculations  and  are  described  in 
Appendix  A: 

(1 )  Geocentric-Equatorial  Inertial  Coordinate  System  of  the  Year  2000  (GEI J2000), 

(2)  Geogr^hic  Coordinate  System  (GEO), 

(3)  Horizon  Coordinate  System  ("South,  East,  and  Zenith"  (SEZ)), 

(4)  Satellite-Fixed  Coordinate  System  (SF), 

(5)  Satellite  Geodetic-Reference  Coordinate  System  (SGR), 

(6)  Geomagnetic  Coordinate  System  (MAG). 

The  convention  used  to  format  this  document  is  illustrated  in  Figure  1.  Definitions  are  followed  by  an 
example,  diagram,  or  equation  only  when  necessary  for  clarification.  AU  italicized  words  are  defined  in 
the  glossary  in  Appendix  B.  The  equations  and  algorithms  for  computing  various  attributes  are  in 
Appendix  C.  Coordinate  transformations  firom  GEI  to  GEO  and  GEO  to  MAG  are  in  Appendix  D. 


RDBMS  TABLE  NAME 


Attribute  Name 

A  short  Description. 

Type:  C.  FORTRAN,  IDL 

Example:  An  example  is  not  given  for  every  attribute. 

Figure:  A  diagram  is  not  given  for  every  attribute. 

Equation:  An  equation  is  not  given  for  every  attribute. 


Figure.  1  -  Document  fonnat 
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2.  SCANS  DATA  TABLE 


2.1  Orbit  File  Name 

An  airay  of  1 1  characters  which  contains  the  name  of  the  orbit  file.  This  name  does  not  include 
device  names  or  subdirectories.  The  first  two  characters  identify  the  database  that  the  file  belongs 
to.  For  RAIDS  this  will  always  be  "ra"  to  show  that  this  file  belongs  to  the  RAIDS  Database.  The 
next  six  characters  contain  the  orbit  number.  The  last  four  characters  will  always  be  ".orb" 
because  this  is  an  orbit  file. 

Type:  unsigned  char[l  1],  CHARACTER*!  1,  String. 

Example:  raOOlOl  .orb  is  the  name  of  the  101“  orbit  file. 


2.2  Beginmng  Record  Byte  Number 

A  number  representing  the  byte  offset  from  the  beginning  of  the  orbit  file  to  the  byte  at  which  the 
first  record  for  the  current  scan  data  begins.  The  first  byte  of  the  file  is  considered  to  be  byte  zero. 
This  follows  the  standard  C  and  IDL  definitions  of  arrays.  The  Beginning  Record  Byte  Number 
actually  locates  the  first  byte  of  the  last  General  Houseteeping  Record  (GHR)^  that  occurred  prior 
to  the  start  of  the  current  scan  data  (see  introduction). 

Type:  long,  INTEGER*4,  Long. 


23  Ending  Record  Byte  Number 

A  number  representing  the  byte  offset  fix)m  the  beginning  of  the  orbit  file  to  the  "end"  of  die  last 
record  for  die  current  scan.  The  first  byte  of  the  file  is  considered  to  be  byte  zero.  If  the  Ending 
Record  Bjde  number  is  less  than  the  Beginning  Record  Byte  Number,  then  the  scan  data  extend 
over  two  orbit  files.  The  byte  offset  would  then  be  fixim  the  beginning  of  the  next  orbit  file.  The 
Ending  Record  Byte  Number  actually  locates  the  last  byte  of  die  first  GHR  that  occurred  after  the 
end  of  the  current  scan  data. 

Type:  long,  INTEGER*4,  Long. 


2.4  Calendar  Year 
The  calendar  year. 

Type:  short,  INTEGER*!,  Integer. 

Example:  1993  is  the  year  1993;  2025  is  the  year  2025. 
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2.5  Day  of  Year 


A  number  from  1  to  366  which  represents  the  number  of  days,  including  the  current  day,  since  the 
beginning  of  the  year.  Note  that  a  month  and  day  can  produce  a  different  Day  of  Year  if  Calendar 
Year  is  a  leap  year. 

Type:  short,  INTEGER*2,  Integer. 

Example:  20  is  January  20;  366  is  December  31  if  Calendar  Year  is  a  leap  year. 


2.6  GMT  Milliseconds 

A  number  from  0  to  86,399,999  which  represents  the  number  of  milliseconds  that  have  elapsed 
since  midnight  Greenwich  Mean  Time  (local  mean  solar  time  on  the  Greenwich  meridian). 

Type:  long,  INTEGER*4,  Long. 

Example:  1,000  is  one  second  after  midnight 


2.7  Scan  Mode  Number 


A  number  equal  to  7  or  a  number  from  0  to  4  which  indicates  the  RAIDS  mode  of  operation  for 
the  present  scan.  The  modes  are  summarized  in  Table  1.  See  the  RAIDS  Experiment  Users 
Manuaf  for  further  information  about  modes  of  operation. 


Type:  unsigned  char,  BYTE,  Byte. 

Example:  0  indicates  that  all  instruments  are  powered  off  and  that  RAIDS  is  waiting  for  the 

next  command. 


Mode  Number 

Characteristic 

Scan  Platform 

Spectrometers 

0 

Power-down  standby 

n/a 

n/a 

1 

Monochromatic  limb  scans 

Continuous  scans 

Fixed  wavelength 

2 

Wavelength  scans  at 
selected  altitudes 

Preset  altitudes 

Scanning 

3 

Fixed  altitude 
wavelength  scans 

Fixed  altitude 

Scanning 

4 

Stare 

Fixed  altitude 

Fixed  wavelength 

7 

Power-up  standby 

n/a 

n/a 

2.8  Instrument  Status 


A  sequence  of  8  toggle  bits  that  form  a  number  from  0  to  255  and  indicates  which  of  the  eight 
RAIDS  instruments  are  turned  on  and  possibly  taking  data.  A  0  in  the  bit's  location  indicates  that 
the  instrument  is  off.  Table  2  shows  which  RAIDS  instrument  corresponds  to  which  bit. 

Type:  unsigned  char,  BYTE,  Byte. 

Example:  255  means  that  all  instruments  are  turned  on;  128  means  that  only  the  EUV 

instrument  is  turned  on. 


Table  2  -  Instrument  Status  Bit  Values 


Bit  Value 

Extreme  Ultraviolet  Spectrograph 

128 

Far  Ultraviolet  Spectrograph 

64 

Middle  Ultraviolet  Spectrometer 

32 

Near  Ultraviolet  Spectrometer 

16 

Near  Infrared  Spectrometer 

8 

5890  A  Photometer 

4 

6300  A  Photometer 

2 

7774  A  Photometer 

1 

2.9  Status  Verified 

A  sequence  of  8  toggle  bits  that  form  a  number  from  0  to  255  which  indicates  that  the  stams  as 
reported  in  Instrument  Status  has  been  verified  by  a  scientist  or  engineer.  A  0  in  the  bit's  location 
indicates  that  the  status  of  that  instrument  has  not  been  verified.  Bits  correspond  to  instruments  as 
described  in  the  table  above. 

Type:  unsigned  char,  BYTE,  Byte. 

Example:  88  means  that  the  FUV,  NUV,  and  NIR  instrument  statuses  have  been 

verified. 


2.10  Image  ID 

This  number  is  generated  by  the  RDBMS  and  points  to  an  image  table. 
Type:  long,  INTEGER’'‘4,  Long. 
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2.11  lonosonde/Radar  Flag 

A  flag  that  indicates  if  any  lonosonde  or  Radar  sites  are  “correlated”  with  the  present  scan.  To  be 
correlated,  the  site  must  be  sufficiently  near  to  the  track  of  the  subtangent  point  during  the  scan. 

Type:  char,  LOGICAL,  Byte. 


2.12  UV  Stars  Flag 

A  flag  that  indicates  if  any  stars  in  the  RAIDS  catalog  are  “correlated”  with  the  present  scan.  The 
stars  in  the  catalog  comprise  all  stars  known  to  emit  strongly  within  the  ultraviolet  wavelength 
range  covered  by  RAIDS.  To  be  correlated  with  the  scan,  the  star’s  right  ascension  (RA)  and 
declination  (DEC)  must  be  sufficiently  near  to  the  RA  and  DEC  ranges  of  the  instrument’s  look 
direction  during  the  scan. 

Type:  char,  LOGICAL,  Byte. 
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3.  EPHEMERIS  AND  MODE  DATA  TABLE 


3.1  GEIJ2000X 

The  X  coordinate  in  centimeters  of  the  satellite's  position  in  the  GEI 12000  system. 

Type:  long,  INTEGER*4,  Long. 

Figure:  GEI  J2000  Coordinate  System,  Appendix  A. 

3.2  GEIJ2000Y 

The  Y  coordinate  in  centimeters  of  the  satellite's  position  in  die  GEI  J2000  system. 

Type:  long,  INTEGER''‘4,  Long. 

Figure:  GEI  J2000  Coordinate  System,  Appendix  A. 

3JGEIJ2000Z 

The  Z  coordinate  in  centimeters  of  the  satellite's  position  in  the  GEI  J2000  system. 

Type:  long,  INTEGER*4,  Long. 

Figure:  GEI  J2000  Coordinate  System,  Appendix  A. 

3.4  Satellite  Altitude 

The  distance  in  centimeters  from  die  subsatelttte  point,  on  the  eardi's  surface,  to  the  satellite. 

Type:  long,  INTEGER’'‘4,  Long. 

Figure:  2 

Equation:  Appendix  C.1 


3.5  Satellite  Distance  fr'om  Earth's  Center 
The  magnitude  in  centimeters  of  the  satellite  vector. 
Type:  long,  INTEGER''‘4,  Long. 
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3.6  Satellite  Geodetic  Latitude 


A  number  from  -9000  to  9000,  in  hundredths  of  a  degree,  which  represents  the  angle  measured 
between  the  equatorial  plane  and  the  downward-projected  normal  to  the  surface  of  the  reference 
ellipsoid  at  the  subsatellite  point.  The  angle  is  measured  from  the  equatorial  plane  with  positive 
angles  being  northward. 

Type:  short,  INTEGER*2,  Integer. 

Example:  4525  is  45.25°  North  Latitude;  -9000  is  the  South  Pole. 

Figiue:  2 

Equation:  Appendix  C.l 


flgnre  2  -  Satellite  altitude  and  geodetic  latitude 
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3.7  Satellite  Longitude 


A  number  from  -17999  to  18000  in  himdredths  of  a  degree  which  represents  the  angle  measured 
positive  eastward  from  the  Greenwich  meridian  to  the  projection  of  fee  satellite  vector  onto  fee 
equatorial  plane. 

Type:  short,  INTEGER*2,  Integer. 

Example:  12135  is  121.35°  East  Longitude;  -500  is  5.00°  West  Longitude. 

Figure:  3 

Equation:  Appendix  C.2 


Figure  3  -  Satellite  longitude  and  local  time.  Note  that  satellite  longitude  is  negative  in  this  figure  because  it  is 
being  measured  westward. 
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3.8  SateUite  Local  Time 


A  number  from  0  to  35999  in  hundredths  of  a  degree  which  represents  the  satellite's  longitude 
relative  to  "midmght".  To  find  this  angle  the  satellite  vector  and  a  vector  with  the  opposite 
direction  as  the  sun  vector  must  both  be  projected  onto  the  equatorial  plane.  The  satellite  local 
time  is  then  measured  eastward  from  the  projection  of  the  sun  vector  to  the  projection  of  the 
satellite  vector. 


Type:  long,  INTEGER*4,  Long. 

Example:  18000  means  that  the  satellite  is  at  1200  hours  local  time  or  noon. 

Figure:  3 

Equation:  Appendix  C.3 


3.9  Satellite  Magnetic  Latitude 

A  number  from  -9000  to  9000  in  hundredths  of  a  degree  which  represents  the  angle  measured 
positive  northward  from  the  magnetic  equatorial  plane  to  the  satellite  vector.  Note  that  this  angle 
is  independent  of  the  shape  of  the  earth. 

Type:  short,  INTEGER*!,  Integer. 

Example:  0  is  the  Magnetic  Equatorial  Plane;  -9000  is  the  Magnetic  South  Pole. 

Figure:  4 

Equation:  Appendix  C.4 


3.10  Satellite  Magnetic  Longitude 

A  number  from  -17999  to  18000  in  hundredths  of  a  degree  which  represents  the  anglp.  measured 
positive  eastward  from  the  +Xmag  of  the  MAG  system  to  the  projection  of  the  satellite  vector 
onto  the  magnetic  equatorial  plane. 

Type: 

Example: 

Figure: 

Equation: 


short,  INTEGER*!,  Integer. 

-5378  is  53.78°  West  Magnetic  Longitude. 
4 

Appendix  C.5 
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Figure  4  -  Satellite  magnetic  ladtutfe  and  longitude 
3.11  Satellite  Solar  Azimuth  Angle 

A  number  from  0  to  35999  in  hundredths  of  a  degree  which  represents  the  angle  measured 
eastward  from  "north''  (-Xu  axis  of  SEZ  system  centered  at  the  satellite)  to  the  projection  of  the 
satellite-sun  vector  onto  the  plane  of  the  horizon  (XhYh  plane). 

Type:  long,  INTEGER’'‘4,  Long. 

Example:  18000  means  that  the  sun  is  directly  south  of  the  satellite. 

Figure:  5 

Equation:  Appendix  C.6 
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3.12  Satellite  Solar  Zenith  Angle 


A  number  from  0  to  18000  in  hundredths  of  a  degree  which  represents  the  angle  from  the  local 
vertical  (Zjj  axis  of  the  SEZ  system  centered  at  the  satellite)  to  the  satellite-sun  vector. 

Type:  short,  INTEGER*2,  Integer. 

Example:  12031  means  the  sun  is  120.31°  from  the  local  vertical,  30.31°  below  the  local 

horizontal. 

Figure:  5 

Equation:  Appendix  C.6 


Figure  5  -  Satellite  solar  azimuth  and  zenith  angles  or  satellite  lunar  azimuth  and 
zenith  angles.  Note  that  vector  s  is  either  the  satellite-sun  vector  or  the 
satellite-moon  vector,  and  the  XhYh  plane  is  parallel  to  the  plane  that  is 
tangent  to  the  Earth's  surface  at  the  subsatellite  point. 


12 


3.13  Satellite  Lunar  Azimuth  Angle 

The  same  as  Satellite  Solar  Azimuth  Angle,  except  that  the  projection  of  the  satellite-moon  vector 
onto  the  plane  of  the  horizon  is  used. 

Type:  long,  INTEGER*4,  Long. 

Figure:  5 

Equation:  Appendix  C.7 


3.14  Satellite  Lunar  Zenith  Angle 

The  same  as  Satellite  Solar  Zenith  Angle,  except  that  the  angle  is  measured  between  the  local 
vertical  and  the  satellite-moon  vector. 

Type:  short,  INTEGER*2,  hiteger. 

Figure:  5 

Equation:  Appendix  C.7 


3.15  Satellite  Pitch 

A  number  between  -17999  and  18000  in  hundredths  of  a  degree  which  represents  the  angle  that 
the  satellite  has  rotated  about  the  +Zg  axis  of  the  SGR  system.  The  right  hand  rule  is  used  to 
define  the  positive  sense  of  this  angle.  Satellite  pitch,  roll,  and  yaw  are  linearly  interpolated  be¬ 
tween  values  given  once  a  second  in  the  Orbital  Parameter  Block  of  the  Circular  File  Record*. 
This  interpolation  is  based  on  the  value  of  GMT  Milliseconds  which  is  defined  above  in  the  Scans 
Data  Table.  The  "once  a  second"  values  of  pitch,  roll,  and  yaw  used  in  the  interpolation  are 
known  to  a  tenth  of  a  degree  in  the  case  of  a  TIROS-N  host 

Type:  short,  INTEGER*2,  Integer. 

Figure:  6 


3.16  Satellite  RoU 

The  same  as  Satellite  Pitch,  except  that  this  angle  represents  the  rotation  about  the  +Ye  axis. 

Type:  short,  INTEGER’''2,  Integer. 

Figure:  6 
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3.17  Satellite  Yaw 


The  same  as  Satellite  Pitch,  except  that  this  angle  represents  the  rotation  about  the  +X(-  axis. 

Type:  short,  INTEGER*2,  Integer. 

Figure:  6 


Xg 

Figure  6  -  Satellite  pitch,  roll,  yaw 

3.18  Satellite  Velocity  X 

The  rate  of  change  in  centimeters  per  second  of  the  satellite's  GEI J2000  X-coordinate. 
Type:  long,  INTEGER*4,  Long 

3.19  Satellite  Velocity  Y 

The  rate  of  change  in  centimeters  per  second  of  the  satellite's  GEI  J2000  Y-coordinate. 
Type:  long,  INTEGER*4,  Long. 
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3.20  Satellite  Velocity  Z 

The  rate  of  change  in  centimeters  per  second  of  the  satellite's  GEI J2000  Z-  coordinate. 

Type:  long,  INTEGER*4,  Long. 

3.21  a  Step  Number 

A  number  between  -150  and  8000  which  represents  the  value  of  the  platform  position  counter. 
Type:  short,  INTEGER*2,  Integer. 


3.22  a 

A  number  nominally  between  1000  and  2650  in  hundredths  of  a  degree  which  represents  the  angle 
that  the  line  of  sight  vector  makes  with  the  plane  of  the  SF  system.  The  angle  is  positive 

when  measured  "earthward"  from  the  satellite  horizontal.  This  value  is  calculated  based  on  the  - 
Step  Number  and  the  platform  calibration  function  determined  before  launch.  Note  that  -  is  only 
one  factor  that  determines  Tangent  Altitude. 

Type:  short,  INTEGER'''2,  Integer. 

Example:  1532  means  that  the  line  of  sight  is  15.32'“  off  the  YbZb  plane  in  the  +Xb 

direction. 


3J13  Tangent  Altitude 

The  distance  in  centimeters  from  die  subtangent  point  to  the  tangent  point. 

Type:  long,  INTEGER*4,  Long. 

Figure:  7 

Equation:  Appendix  C.8 


Figure  7  -  Tangent  altitude 


3^  Tangent  Geodetic  Latitude 

The  same  as  Satellite  Geodetic  Latitude,  excq)t  that  the  angle  is  measured  from  the  equatorial 
plane  to  the  downward-projected  normal  of  the  reference  ellipsoid  at  the  subtangent  point. 

Type:  short,  INTEGER*2,  Integer. 

Equation:  Appendix  C.9. 


3.2S  Tangent  Longitude 

The  same  as  Satellite  Longitude,  except  that  die  projection  of  the  tangent  vector  onto  the 
equatorial  plane  is  used. 

Type:  short,  INTEGER*2,  Integer. 

Equation:  Appendix  C.10. 
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3.26  Tangent  Local  Time 


The  same  as  Satellite  Local  Time,  except  that  the  angle  is  measured  to  the  projection  of  the 
tangent  vector  onto  the  equatorial  plane. 

Type:  long,  INTEGER*4,  Long. 

Equation:  Appendix  C.  1 1 . 


3.27  Tangent  Magnetic  Latitude 

The  same  as  Satellite  Magnetic  Latitude,  except  that  the  angle  is  measured  between  the  magnetic 
equatorial  plane  and  the  tangent  vector. 

Type:  short,  INTEGER*2,  Integer. 

Equation:  Appendix  C.12. 


3.28  Tangent  Magnetic  Longitude 

The  same  as  Satellite  Magnetic  Longitude,  except  that  the  projection  of  the  tangent  vector  onto  the 
magnetic  equatorial  plane  is  used. 

Type:  short,  INTEGER*2,  Integer. 

Equation:  Appendix  C.13 


3.29  Tangent  Solar  Azimuth  Angle 

The  same  as  Satellite  Solar  Azimuth  Angle,  except  that  the  SEZ  system  is  centered  at  the  tangent 
point  and  the  projection  of  tangent-sun  vector  onto  the  plane  of  the  horizon  is  used. 

Type:  long,  INTEGER*4,  Long. 


330  Tangent  Solar  Zenith  Angle 

The  same  as  Satellite  Solar  Zenith  Angle,  except  that  the  SEZ  system  is  centered  at  the  tangent 
point  and  die  angle  is  measured  between  the  local  vertical  and  die  tangent-sun  vector. 

Type:  short,  INTEGER*2,  Integer. 


331  Tangent  Lunar  Azimuth  Angle 

The  same  as  Tangent  Solar  Azimuth  Angle,  except  that  the  projection  of  the  tangent-moon  vector 
onto  the  plane  of  the  horizon  is  used. 

Type:  long,  INTEGER*4,  Long. 
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3.32  Tangent  Lunar  Zenith  Angle 

The  same  as  Tangent  Solar  Zenith  Angle,  except  that  the  angle  is  measured  between  the  local 
vertical  and  the  tangent-moon  vector. 

Type:  short,  E'JTEGER*2,  Integer. 


3J3  Line  of  Sight  Azimuth  Angle 

The  same  as  Satellite  Solar  Azimuth  Angle,  except  that  the  projection  of  the  line  of  sight  vector  onto 
the  plane  of  the  horizon  is  used.  The  SEZ  system  is  centered  at  the  satellite. 

Type:  long,  INTEGER*4,  Long. 


334  Line  of  Sight  Zenith  Angle 

The  same  as  Satellite  Solar  Zenith  Angle,  except  that  the  angle  is  measured  between  the  local  vertical 
and  the  line  of  sight  vector.  The  SEZ  system  is  centered  at  the  satellite. 

Type:  short,  INTEGER*2,  Integer. 


335  Solar  Right  Ascension 

A  number  from  0  to  235959  in  hours,  minutes,  and  seconds  (hhmmss)  which  represents  the  angle 
measured  positive  eastward  (+Y  direction)  from  the  +X  axis  of  the  GEI J2000  system  to  the  projection 
of  the  sun  vector  onto  the  equatorial  plane. 

Type:  long,  INTEGER*4,  Long. 

Example:  123447  is  12  hours,  34  minutes,  and  47  seconds. 

Figure:  8 

Equation:  Appendix  C.  1 4. 


336  Solar  Declination 

A  number  from  -900000  to  900000  in  degrees,  minutes,  and  seconds  (°° ' ' ''  ")  which  represents  the 
angle  measured  positive  northward  (+Z  direction  of  the  GEI  J2000  system)  from  the  equatorial  plane 
to  the  sun  vector. 

Type:  long,  INTEGER*4,  Long. 

Example:  -872358  is  -87  degrees,  23  minutes,  and  58  seconds. 

Figure:  8 

Equation:  Appendix  C.14. 
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3J7  Lunar  Right  Ascension 

The  same  as  Solar  Right  Ascension,  except  that  the  projection  of  the  moon  vector  onto  the  equatorial 
plane  is  used. 

Type:  long,  INTEGER*4,  Long. 

Equation:  Appendix  C.  15. 


338  Lunar  Declination 

The  same  as  Solar  Declination,  except  that  the  angle  is  measured  between  the  equatorial  plane  and 
the  moon  vector. 

Type:  long,  INTEGER*4,  Long. 

Equation:  Appendix  C.15 


Figure  8  -  Solar  right  ascension  and  declination 
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3.39  Satellite  Right  Ascension 


The  same  as  Solar  Right  Ascension,  except  that  the  projection  of  the  satellite  vector  onto  the 
equatorial  plane  is  used. 

Type:  long,  INTEGER*4,  Long. 

3.40  Satellite  Declination 

The  same  as  Solar  Declination,  except  that  the  angle  is  measured  between  the  equatorial  plane  and 
the  satellite  vector. 

Type:  long,  INTEGER*4,  Long. 


3.41  Line  of  Sight  Right  Ascension 

The  same  as  Solar  Right  Ascension,  except  that  a  vector  starting  at  the  GEI J2000  origin  with  the 
same  direction  as  the  Une  of  sight  vector  is  projected  onto  the  equatorial  plane. 

Type:  long,  INTEGER*4,  Long. 

3.42  Line  of  Sight  Declination 

The  same  as  Solar  Declination,  except  that  the  angle  is  measured  between  the  equatorial  plane  and 
a  vector  that  has  the  same  direction  as  the  line  of  sight  vector  and  that  starts  at  the  GEI  J2000 
origin. 

Type:  long,  INTEGER*4,  Long. 


3.43  Line  of  Sight  to  Sun  Angle 

A  number  between  0  and  18000  in  hundredths  of  a  degree  which  represents  the  angle  measured 
between  the  Une  of  sight  vector  and  the  satelUte-sun  vector. 

Type:  short,  INTEGER*2,  Integer. 

Figure:  9 

3.44  Line  of  Sight  to  Moon  Angle 

The  same  as  Line  of  Sight  to  Sun  Angle,  except  that  the  angle  is  measured  between  the  line  of 
sight  vector  and  the  satellite-moon  vector. 

Type:  short,  INTEGER*2,  Integer. 


20 


Figure  9  •  Line  of  sight  to  sun  angle 


3.45  MUV  X.  Step  Number 

A  number  between  -1  and  2000  which  indicates  the  grating  position  step  coimt  of  the  MUV 
spectrometer.  A  value  of  -1  means  that  the  spectrometer  grating  is  scanning. 

Type:  short,  INTEGER*2,  Integer. 

3.46  NUV  k  Step  Number 

The  same  as  MUV  k  Step  Number,  except  that  this  is  for  the  NUV  spectrometer. 

Type:  short,  INTEGER*2,  hiteger. 

3.47  NIR  k  Step  Number 

The  same  as  MUV  A  Step  Nmnber,  except  that  this  is  for  the  NIR  spectrometer. 

Type:  short,  INTEGER*2,  Integer. 
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3.48  MUVA 


A  number  in  angstroms  which  indicates  the  MUV  spectrometer  wavelength  setting.  This  value  is 
calculated  based  on  the  MUV  X  Step  Number  and  the  ciurent  wavelength  calibration  function.  A 
value  of  -1  means  that  the  spectrometer  grating  is  scanning. 

Type:  short,  INTEGER*2,  Integer. 

Example:  2450  means  that  the  spectrometer  is  looking  at  2450  angstroms. 


3.49  NUVA 

The  same  as  MUV  X,  except  that  this  is  for  the  NUV  spectrometer. 
Type:  short,  INTEGER*2,  Integer. 


3.50  NIR  X 

The  same  as  MUV  X,  except  that  this  is  for  the  NIR  spectrometer. 
Type:  short,  INTEGER*2,  Integer. 


3.51  EUV  Grating  Position 

A  toggle  that  indicates  the  position  of  the  Extreme  Ultraviolet  Spectrograph's  grating.  Zero 
indicates  that  the  grating  is  positioned  for  the  shorter  wavelength  setting;  any  other  value  indicates 
the  longer  wavelength  setting. 

Type:  unsigned  char,  BYTE,  Byte. 


3.52  Day/Night 

A  toggle  that  indicates  the  status  of  the  RAIDS  day/night  clock.  Zero  indicates  day;  any  other 
value  indicates  night.  If  RAIDS  has  been  issued  a  sequence  command,  then  the  day/night  toggle 
provides  information  as  to  which  part  of  the  sequence  RAIDS  is  currently  executing.  See  the 
RAIDS  Experiment  Users  Monuaf  for  ftirther  information  about  sequencer  operation. 

Type:  unsigned  char,  BYTE,  Byte. 
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Appendix  A 

COORDINATE  SYSTEMS®^*® 


A.1  Geocentric-Equatorial  Inertial  Coordinate  System  of  the  Year  2000 


The  Geo^ntric-Equatorial  Inertial  Coordinate  System  of  the  Year  2000  (GEI J2000)  is  an  inertial  system 
whose  origin  is  located  at  the  center  of  the  earth.  The  fundamental  plane  is  fliat  of  the  earth's  geographic 
equator.  The  positive  X  axis  points  in  the  direction  of  the  ascending  node  of  the  ecliptic  (mean  plane  of 
the  earth's  orbit  around  the  Sun)  on  the  earth's  geographic  equator  at  noon  barycentric  dynamical  time  on 
January  1, 2000  (Julian  day  number  2,451445.0).  The  positive  Z  axis  points  in  the  direction  of  the  earth's 
north  pole,  and  the  Y  axis  completes  the  right-handed  orthogonal  set. 
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A.2  Geographic  Coordinate  System 


\ 

X 


The  geogr^hic  coordinate  system  (GEO)  is  a  rotating  system  whose  origin  is  located  at  the  center 
of  the  earth.  The  fundamental  plane  is  that  of  the  earth's  geographic  equator.  The  positive  X  axis  is  fixed 
with  the  rotation  of  the  earth  so  that  it  passes  through  the  Greenwich  meridian.  The  positive  Z  axis  points 
in  the  direction  of  the  earth's  north  pole,  and  the  Y  axis  completes  the  right-handed  orthogonal  set. 
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A3  Horizon  Coordinate  System 


The  Horizcm  Coordinate  System  (SE^  is  commonly  referred  to  as  a  topocentric-horizon  system. 
This  means  that  the  origin  of  the  system  is  located  at  the  observer,  i.e.  RAIDS.  The  local  vertical,  located 
on  the  near  side  of  the  reference  ellipsoid  and  passing  through  die  observer,  defines  the  Zg  axis;  the 
positive  sense  is  away  from  the  earth.  The  positive  Xh  axis  points  due  south,  while  the  positive  Yh  axis 
points  due  east.  The  directions  for  east  and  south  are  interpreted  as  if  the  origin  were  on  the  earth's  surface. 
The  SEZ  system  is  not  an  inertial  reference  fiame. 
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A.4 


Satellite-Fixed  Coordinate  System 


Velocity 
^Vector 


The  Satellite-Fixed  Coordinate  System  (SF)  is  a  right-handed,  body  fixed-coordinate  system.  The 
SF  coordinate  system  is  defined  on  page  8-1  in  reference  3  as  follows: 

The  Zg  (pitch)  axis  is  normal  to  the  separation  plane  and  passes  through  the  center  of  the 
separation  ring;  the  positive  sense  of  the  Zg  axis  is  nominally  along  the  positive  orbit  normal.  The 
Xb  (yaw)  axis  is  normal  to  the  earth-facing  side  of  the  spacecraft  with  the  positive  sense  being 
directed  towards  the  earth.  The  Yg  (roll)  axis  completes  the  right-handed  orthogonal  triad;  the 
positive  sense  of  the  roll  axis  is  nominally  opposite  to  the  instantaneous  velocity  vector  of  the 
satellite.  The  origin  of  the  reference  coordinate  system  is  located  on  the  plane  of  separation 
between  the  spacecraft  and  its  booster. 

Note  that  the  phrase  "reference  coordinate  system"  refers  to  the  SF  system.  Given  that  RAIDS  will 
probably  fly  on  a  satellite  other  than  a  NOAA  TIROS-N,  some  of  the  terminology  of  this  description  might 
change.  However,  the  definition  shown  by  the  above  figure  will  likely  remain  the  same. 
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A.5  Satellite  Geodetic-Reference  Coordinate  System 


The  Satellite  Geodetic-Reference  Coordinate  System  (SGR)  is  also  defined  in  reference  3  as 

follows: 

The  local  geodetic  reference  represents  the  nominal  required  directions  for  the  spacecraft 
coordinate  reference.  The  Xq  axis  is  a  line  through  the  origin  of  the  spacecraft  coordinate 
reference  and  normal  to  the  reference  geoid  on  the  near  side.  Xq  is  positive  from  the  spacecraft 
toward  the  earth.  The  axis  is  a  line  normal  to  the  Xq  axis  and  to  the  instantaneous  inertial 
velocity  vector.  The  positive  Zq  axis  lies  on  the  sunlit  side  of  die  plane  formed  by  the  velocity 
vector  and  the  Xq  axis.  The  Yq  axis  completes  the  right-handed  orthogonal  set 

Note  that  'local  geodetic  reference"  refers  to  the  SGR  system  and  that  "spacecraft  coordinate  reference" 
refers  to  the  SF  system. 


A.6  Geomagnetic  Coordinate  System 


The  Greomagnedc  Coordinate  System  (MAG)  is  a  non-inertial  system  that  rotates  with 
the  earth.  This  paper  assumes  a  Centered  Dipole  Geomagnetic  coordinate  system.  Its  origin  is  located  at 
the  geographic  center  of  the  earth  with  the  2^g  axis  running  along  the  earth's  magnetic  dipole.  The 
positive  sense  of  the  Z^ag  axis  is  towards  the  north  pole  from  the  center  of  the  earth.  The  positive  Y^ag 
axis  is  defined  as  the  cross  product,  x  Zmag,  of  the  geographic  north  pole  (+Z  GEI J2000  axis)  and  the 
Zj4ag  axis.  Finally,  the  positive  axis  is  defined  as  x  Z^ag.  Note  that  the  Xi^qYmag  plane 
defines  the  magnetic  equator. 
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Appendix  B 
GLOSSARY 


Circular  File 

Line  of  Sight  Vector 

Magnetic  Equatorial 
Plane 

Moon  Vector 
Reference  Ellipsoid 

Satellite  Vector 

Satellite>Moon  Vector 

Satellite*Sun  Vector 

Subsatellite  Point 

Subsatellite  Vector 
Subtangent  Point 

Subtangent  Vector 


NOAA  file  containing  the  sateUite  data  and  RAIDS  instrument  data  for 
ingest  into  RAIDS  Data  Processing  Facility  (DPF).  The  data  will  be 
transferred  from  the  NOAA  Circular  File  to  the  RAIDS  DPF  over 
Internet,  one  pass  (orbit)  at  a  time. 

The  vector  pointing  from  RAIDS  in  the  direction  of  the  field  of  view  of 
the  instruments  mounted  on  the  scanning  platform. 

The  plane  containing  the  X^ao  and  Y^ao  axes  of  the 
Geomagnetic  Coordinate  System  (MAG).  See  appendix  A  for 
information  about  coordinate  systems.  Note  that  this  plane  passes 
through  the  center  of  the  earth. 

The  vector  from  the  center  of  the  earth  to  the  center  of  the  moon. 

A  model  that  t^proximates  the  surface  of  the  earth,  commonly  referred  to 
as  "mean  sea  level"  or  "geoid".  Values  used  for  this  model  come  from 
The  Astronomical  Almanac,  page  K5^. 

The  vector  from  the  center  of  the  earth  to  the  satellite's  position.  This 
vector  has  components  equal  to  the  GEI J2000  X,  Y,  and  Z  coordinates. 

The  vector  from  the  satellite's  position  (origin  of  Satellite-Fixed 
coordinate  system,  Appendix  A4)  to  the  center  of  the  moon. 

The  vector  from  the  satellite's  position  (origin  of  SateUite-Fixed 
coordinate  system,  Appendix  A.4)  to  the  center  of  the  sun. 

The  point  on  the  r^erence  ellipsoid  where  the  outward  normal  to  the 
reference  ellipsoid  passes  through  the  sateUite's  position. 

The  vector  from  the  center  of  the  earth  to  the  subsatellite  point.  See 
Figure  2,  page  8. 

The  point  on  the  reference  ellipsoid  where  the  local  vertical  is  perpendic¬ 
ular  to  the  line  of  sight  vector  and  for  which  the  distance  to  the  line  of 
sight  vector  is  a  minimum. 

The  vector  from  the  center  of  the  earth  to  die  subtangent  point. 
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Sun  Vector 


The  vector  from  the  center  of  the  earth  to  the  center  of  the  sun. 


Tangent  Point 


Tangent-Moon  Vector 
Tangent-Sun  Vector 
Tangent  Vector 


The  point  along  the  line  of  sight  vector  where  the  distance  between  the 
line  of  sight  vector  and  the  reference  ellipsoid  is  a  minimum.  The  local 
vertical  at  the  sub  tangent  point  passes  through  the  tangent  point. 

The  vector  from  the  tangent  point  to  the  center  of  the  moon. 

The  vector  from  the  tangent  point  to  the  center  of  the  sum. 

The  vector  from  the  center  of  the  earth  to  the  tangent  point. 
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Appendix  C 
EQUATIONS 

The  equations  and  algorithms  below  are  central  to  computing  various  attributes  described  in  Sections  2  and 
3.  Some  aspects  of  the  notation  and  the  mathematical  functions  require  further  specification.  The  notation 
for  unit  vectors  in  any  coordinate  system  is  y,  and  £.  If  a  given  attribute  (e.g.  Satellite  altitude)  is 
discussed  in  the  context  of  a  particular  coordinate  system,  then  the  reader  should  assume  that  i,  jp, 
and  t  are  the  unit  vectors  in  that  coordinate  system.  The  definitions  of  and  t  can  change  from 
attribute  to  attribute.  In  aU  cases,  the  equations  below  were  written  to  directly  correspond  with  the  actual  C 
code. 

Several  equations  listed  below  implement  the  arc  tangent  (tan'')  function.  The  arc  tangent  (tan'')  is 
computed  with  the  C  function  atan2.  This  function  returns  a  value  between  -it  and  it  radians.  It  computes 
atan2(r]/ro)  using  the  signs  of  r^  and  Iq  to  determine  the  quadrant  of  the  arc  tangent.  Atan2  is  well  defined 
when  either  r,  or  rj  is  zero,  but  may  flag  a  domain  error  if  both  arguments  are  zero'^ 


C.1  Satellite  Altitude  and  Geodetic  Latitude* 

The  following  items  are  required  for  the  calculation: 

^  ''r  S  *  t  (1) 


(2) 


where  Rgq  is  a  three  dimensional  GEI  vector  in  centimeters  from  the  center  of  the  Earth  to  the  satellite. 
Define 


1-(2/-/)cos24)' 


6  =  tan"* 


l)/^J 


.  -90'’:s6i90' 


(3) 


(4) 


where: 

a^  -  the  eqtiatorial  radius  of  the  Earth -637814000.0  cm. 
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f  -  the  flattening  factor  of  the  Earth  «  1 .0/298.257. 

4),'  -  the  geocentric  latitude  of  the  subsatellite  point. 

At  this  step,  (j);  is  initially  assumed  to  equal  6.  The  geodetic  latitude  (<t)J  and  altitude  (A),  respectively  are 
as  follows: 


/ 

(J)^  =  tan'^ 

k 


1 

(1-/)' 


-90°^  4)^  ^90° 


(5) 


A  =  -  r^cos(4)^-4)')  • 


(6) 


Given  the  above  quantities,  define  updated  quantities  via  equations  (7)  and  (8): 


A4): 


sm 


-sm(4) -4)') 


4)'  =  6  -  A4)'  . 


(7) 

(8) 


Equations  (3)  -  (8)  are  repeated  until  the  new  4),'  is  within  1  x  10  ‘°  of  the  previous  value.  When  this 
occurs  4)s  is  the  calculated  geodetic  latitude  in  radians  and  A  is  the  altitude  in  centimeters.  The  geodetic 
latitude  is  then  converted  to  hundredths  of  a  degree  by  multiplying  by  ISOOO/tc  and  rounding  to  the  nearest 
whole  number.  If  this  loop  is  performed  1000  times  it  will  cause  the  function  to  stop  and  issue  an  error 
message. 


Satellite  Longitude 

The  following  vector  is  required  for  the  calculation: 


'^GEO 


=  r.-x  +  rry  +  r,-f 


where  R^eo  is  a  three  dimensional  GEO  vector  from  the  center  of  the  Earth  to  the  satellite. 
Calculate  the  satellite  longitude  (satjlong)  as  follows: 


sot-long  = 


0.0 

for 

'■o  =  '■l  =  ' 

um-f  :■  ]  •  M 

otherwise. 

V  ''oy 


It 
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Longitude  is  converted  to  hundredths  of  a  degree  by  multiplying  by  100  and  rounding  to  the  nearest  whole 
number. 

Satellite  Local  Time 

The  following  items  are  required  for  the  calculation: 

ra^  -  right  ascension  of  the  satellite  in  degrees. 

right  ascension  of  the  sun  in  degrees. 

Calculate  the  satellite  local  time  as  follows: 

=  Kor  -'■«««)-  180-0 

If  the  computed  satellite  local  time  is  greater  than  or  equal  to  360.0,  subtract  360.0  to  arrive  at  the 

final  result.  If  satellite  local  time  is  less  than  0.0 ,  add  360.0.  Local  time  is  converted  to  hundredths  of  a 
degree  by  multiplying  by  100  and  rounding  to  the  nearest  whole  number. 


C.4  Satellite  Magnetic  Latitude 

The  following  vector  is  required  for  the  calculation: 

^MAG  =  +  r^'t 


where  Rmag  ts  ^  three-dimensional  MAG  vector  from  the  center  of  the  earth  to  the  satellite. 
Calculate  the  satellite  magnetic  latitude  {satjat^  as  follows: 


sat-lat^^  =  i 

mag 


0.0 


for 


tan 


-1 


.  180 


''o  =  =  '-2  =  0 

otherwise 


Latitude  is  converted  to  hundredths  of  a  degree  by  multiplying  by  100  and  rounding  to  the  nearest  whole 
number. 


C.5  Satellite  Magnetic  Longitude 

The  following  vector  is  required  for  the  calculation: 

^UAG  = 

where  Rm^g  is  a  three-dimensional  MAG  vector  from  the  center  of  the  earth  to  the  satellite. 
Calculate  the  satellite  magnetic  longitude  (satjong^)  as  follows: 


35 


sot-long 


0.0 


180 


for  =  r^  =  0 

otherwise. 


mag 


tan-1 


V 


'o; 


It 


Longitude  is  converted  to  hundredths  of  a  degree  by  multiplying  by  100  and  rounding  to  the  nearest  whole 
number. 


C.6  Satellite  Solar  Azimuth  and  2^nith 
The  following  vector  is  required  for  the  calculation: 


where  Rsez(O)  is  the  SEZ  satellite-sun  vector. 
Calculate  the  solar  azimuth  (JSA)  as  follows: 


SA  = 


0.0 


for 


tan 


-1 


-r. 


0/ 


180 

7C 


To  =  Tj  =  0 

Otherwise. 


If  iSA  is  less  than  0.0,  add  360.0  to  obtain  the  final  result. 
Calculate  the  solar  zenith  (SZ)  as  foDows: 


90.0 

/-T - 

SZ  = 

tan"* 

/  2  2 
h  ^  '■i 

180 

\  2  y 

tu 

for 


'■o  =  ^1  =  ^2  =  0 

otherwise. 


If  SZ  is  less  than  0.0,  add  180.0. 

Solar  azimuth  and  zenith  are  converted  to  himdredths  of  a  degree  by  multiplying  by  100  and  rounding  to 
the  nearest  whole  number. 


C.7  Satellite  Lunar  Azimuth  and  Zenith 
The  following  vector  is  required  for  the  calculation: 

where  Rctt(C)  is  the  SEZ  satellite-moon  vector. 
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Calculate  the  satellite  lunar  azimuth  (LA)  as  follows: 

0.0 

I 

LA  = 


for 


tan' 


t  \ 


180 

1C 


'■o  =  '•■  =  0 


Otherwise. 


If  LA  is  less  than  0.0,  add  360.0  to  obtain  the  final  result. 
Calculate  the  satellite  lunar  zenith  (LZ)  as  follows: 


LZ  = 


90.0 

tan  ‘  -i- 


for 


h  +  ^1 


180 

It 


'■o  =  '■i  =  ''z  =  0 


otherwise. 


If  LZis  less  than  0.0,  add  180.0. 

Lunar  azimuth  and  zenith  are  converted  to  hundredths  of  a  degree  by  multiplying  by  100  and  rounding  to 
the  nearest  whole  number. 


C.8  Tangent  Altitude 

The  following  vectors  are  required  for  the  calculation: 

^GET  “  *  r^'y  +  (1) 

^GB  =  *  nil' t  (2) 


where: 

Rgh  is  a  three  dimensional  GEI  vector  in  centimeters  from  the  center  of  the  Earth  to  the  satellite, 
Mgq  is  the  line  of  sight  vector  in  centimeters. 

Define  unit  vectors: 


f  = 


(3) 


Using  the  definition  of  scalar  product,  the  supplement,  6,  of  the  angle  between  the  vectors  r  and  m  can  be 
found: 
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-cos  6  =  rh'f  . 


(4) 


The  tangent  altitude  unit  vector,  d,  can  now  be  calculated: 


a  = 


f 

sin  6 


+ 


th 

tan  6 


(5) 


The  next  step  is  to  calculate  the  subtangent  point.  The  calculation  yields  two  candidates:  one  for  each 
"side"  of  the  ellipsoid  relative  to  the  line  of  sight.  There  are  three  different  cases  to  consider  when 
calculating  the  subtangent  point  depending  on  which,  if  any,  components  of  the  tangent  altitude  unit  vector 
d  =  (oq,  fli,  flj)  are  zero.  The  last  component  needs  to  be  nonzero. 

Case  1 :  If  Oj  =  0  then  skip  to  equation  (6). 

Case  2:  If  Oj  >  0  and  Oq  “  0;  rotate  array  around  so  last  element  is  nonzero: 

a'o  “  Oi;  0^1  -  02,  t^2  -  Oo-  Perform  equations  (6)  through  (16)  using  (f  in  place  of  a.  After 
equation  (16)  return  here. 

The  calculated  subtangent  point  is  in  rotated  variables,  let  s'o  -  So;  s',  -s,;  -  Sj.  Rotate  the 

array  back  around  into  the  proper  order  Sq  -  s'j;  s,  -  s'q;  Sj  -  s',.  Return  to  equation  (17) 

Case  3:  If  Oj  “  0  '  0: 

Rotate  the  altitude  vector  array  so  the  last  element  is  non  zero; 

“  flol  Perform  equations  (6)  through  (16)  replacing  a  with  After 

equation  (16)  return  here. 

The  calculated  subtangent  point  is  in  rotated  variables,  let  s'o  ■  Sq;  s',  -s,;  -  Sj.  Rotate  the 

array  back  around  into  the  proper  order  Sq  -  s',;  s,  -  s'j;  Sj  -  s'o .  Return  to  equation  (17). 

Calculate  the  altitude  component  ratios,  tIq,  and  T),2: 

%  =  ^  ;  <1.2  =  ^  (« 


Calculate  the  first  subtangent  point  and  the  distance  d,  to  the  satellite:  first  define  the  vector  b,  for 
which  bo  -  b,  -  earth  equatorial  radius  -  637814000.0  cm,  and  bj  -  earth  polar  radius  -  bo(1.0  -  f),  where  f 
is  the  flattening  factor  -  1.0/298.257.  Then  we  have 


*0  2  bl  2  , 

b!  bl 


(7) 
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(8) 


„(!)  ^0 
=  Pi  —,^01 


(9) 


Calculate  the  second  subtangent  point  p*  and  the  distance  to  the  satellite  dj: 

n<2)  _  „(1) 

Pq  ~  ~Po 


pr  =  -p1" 


n<2)  _  _(!) 

Pi  ~  ~P2 


The  shorter  of  the  two  distances  is  the  subtangent  point  we  want: 

if  <  d^;  Sq  Jj  =  pl^\  = 

f  Sq  =po®^;  ^1  =P?\  S^  = 


(10) 

(11) 

(12) 

(13) 

(14) 

(15) 

(16) 


The  subtangent  point  vector,  S,  is  now  known: 

S  =  Sq-  i  *  s^•  S  +  S2-  £ 


(17) 
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Using  the  scalar  product  and  some  manipulation,  the  magnitude  of  the  tangent  altitude  vector,  a,  can  be 
found. 


r; 


s.  - 


m 


(18) 


^  ^m-COS(9) 

sin^(0) 


(19) 


a  =  (R  -  R)  '  sm(0) 


(20) 


C.9  Tangent  Geodetic  Latitude 

The  following  GEI  vector  is  required  for  the  calculation: 

d  =  +  tZj-y  +  Oj-z 

where  d  is  the  tangent  altitude  unit  vector  that  was  found  in  the  Tangent  Altitude  Section,  C.8. 
Calculate  the  tangent  geodetic  latitude  (tan_lat)  as  follows: 


0.0 


tan 


-1 


[f 


0  *  ^\  ) 


180 

71 


for  Aq  =  flj  =  flj  =  0 
otherwise 


Tangent  geodetic  latitude  is  converted  to  hundredths  of  a  degree  by  multiplying  by  100  and  rounding  to  the 
nearest  whole  number. 


C.IO  Tangent  Longitude 

The  following  vector  is  required  for  the  calculation: 

^(XO~  ^2*^ 


where  Tc^o  is  the  GEO  tangent  point  position  vector  in  centimeters.  This  vector  is  found  by  first  aHWino 
the  subtangent  point  vector,  s,  and  the  tangent  altitude  vector,  c,  found  in  the  Tangent  Altitude  Section, 
C.8.  The  resulting  vector  is  the  tangent  point  in  the  GEI  system.  That  must  then  be  transformed  into  the 
GEO  system  by  the  method  in  Appendix  D. 
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Calculate  the  tangent  longitude  {tanjong)  as  follows: 


0.0 


for  fj,  =  =  0 


tan-long  = 


/ 

tan"* 


\ 


180 

7C 


Otherwise 


Longitude  is  converted  to  hundredths  of  a  degree  by  multiplying  by  100  and  rounding  to  the  nearest  whole 
number. 


C.11  Tangent  Local  Time 

The  following  items  are  required  for  the  calculation: 

ra^  -  light  ascension  of  the  tangent  point  in  degrees. 

-  right  ascension  of  the  sun  in  degrees. 

Calculate  the  tangent  local  time  ( as  follows: 

'  180-0 

If  the  computed  tangent  local  time  is  greater  than  or  equal  to  360.0,  subtract  360.0  to  obtain  the  final 

result  If  tangent  local  time  is  less  than  0.0  add  360.0.  Tangent  local  time  is  converted  to  hundredths  of  a 
degree  by  multiplying  by  100  and  roimding  to  the  nearest  whole  number. 


C.12  Tangent  Magnetic  Latitude 

The  following  vector  is  required  for  the  calculation: 

^MAG~ 


where  T^ac  is  the  MAG  tangent  point  position  vector  in  centimeters.  This  vector  is  found  by  adding  the 
subtangent  point  vector  and  the  tangent  altitude  vector  found  in  the  Tangent  Altitude  Section,  C.8.  The 
resulting  vector  is  the  tangent  point  in  the  GEI  system.  That  must  then  be  transformed  into  the  MAG 
system. 

Calculate  the  tangent  magnetic  latitude  as  follows: 


0.0 


tan 


-1 


180 

It 


for  to  =  =  ^2  =  0 

otherwise 


Latitude  is  converted  to  himdredths  of  a  degree  by  multiplying  by  100  and  rounding  to  the  nearest  whole 
number. 
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C.13  Tangent  Magnetic  Longitude 

The  following  vector  is  required  for  the  calculation: 


'^MAG 


= 


where  Tm^g  is  the  MAG  tangent  point  position  vector  in  centimeters.  This  vector  is  found  by  first  adding 
the  subtangent  point  vector  and  the  tangent  altitude  vector  found  in  the  Tangent  Altitude  Section,  C.8. 
The  resulting  vector  is  the  tangent  point  in  the  GEI  system.  That  must  then  be  transformed  into  the  MAG 
system  by  the  method  in  Appendix  D. 


Calculate  the  tangent  magnetic  longitude  itanjong^c)  as  follows: 


0.0 


I  tan'* 


180 

It 


for  ^0  =  =  0 

otherwise 


Longitude  is  converted  to  hundredths  of  a  degree  by  multiplying  by  100  and  rounding  to  the  nearest  whole 
number. 


C.14  Solar  Right  Ascension  and  Declination 
The  following  vector  is  required  for  the  calculation: 

+  r,-y'  +  Tj-  f 


where  Rgii(O)  is  a  three  dimensional  GEI  vector  in  centimeters  firom  the  center  of  the  earth  to  the  center  of 
the  sun.  This  vector  is  calculated  using  the  Jet  Propulsion  Laboratory  (JPL)  Export  Ephemeris  Package, 


Calculate  the  right  ascension  (ra)  as  follows: 

0.0 


ra  =  i 


tan' 


(  \ 


\  ''o; 


180 

7t 


for  Tq  =  r,  =  0 
otherwise 


If  ra  is  less  than  0.0,  add  360.0. 
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Calculate  the  declination  (dec)  as  follows: 


dec 


0.0 


tan 


-1 


U‘ 


''o  +  '•i 


180 

71 


for  Tq  =  Tj  =  r2  =  0 
otherwise 


Solar  right  ascension  and  declination  are  then  converted  to  degrees,  minutes,  seconds  format  from  the 
decimal  degrees  calculated. 


C.15  Lunar  Right  Ascension  and  Declination 
The  following  vector  is  required  for  the  calculation: 

+  '•l-y  +  ^2*  ^ 


where  RqqC  C)  is  a  three  dimensional  GEI  vector  in  centimeters  from  the  center  of  the  earth  to  the  center  of 
the  moon.  This  vector  is  calculated  using  the  Jet  Propulsion  Laboratory  (JPL)  Export  Ephemeris  Package. 

Calculate  the  right  ascension  (ra)  as  follows: 


ra  =  i 


0.0 


tan* 


/  \ 


V 


180 

1C 


If  ra  is  less  than  0.0,  add  360.0. 
Calculate  the  declination  (dec)  as  follows: 


for  Tq  =  r,  =  0 
otherwise 


dec 


0.0 


tan 


-1 


180 

7t 


for  To  =  r,  =  =  0 

otherwise 


Lunar  right  ascension  and  declination  are  then  converted  to  degrees,  minutes,  seconds  format  from  the 
decimal  degrees  calculated. 
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Appendix  D 

COORDINATE  TRANSFORMATIONS’***  ” 

D.1  GEI  To  GEO 

Let  the  angle  from  the  X-axis  of  GEI  to  the  Greenwich  meridian  measured  eastward  in  the  Earth's  equator 
be  0.  The  angle  0  is  referred  to  as  the  Greenwich  Mean  Sidereal  Time. 

The  transfer  from  GEO  to  GEI  is 

cos  0  sin  0  0 
-sin  0  cos  0  0 
0  0  1 

or 

^GEO  =  Kcos®)  %  *  (sin0)  -ly  i  +  [(-sin0)  -V,  +  (cos0)  -Vy]  S  +  [VJ  i 

and  the  inverse  transformation  is 

COS  0  -sin  0  0 
sin  0  cos  0  0 
0  0  1 

or 

^ge7  =  KCOS0)  +  (sin0)  -ty  x  +  [(-sin0)  -C^  +  (cos0)  -ty  y  +  [t/J  z 


GEO  To  MAG 


Let  6  be  the  GEO  colatitude  and  (J)  the  GEO  longitude  of  the  centered  dipole  axis  in  its  extension  through 
the  northern  hemisphere  of  the  earth.  The  transformation  from  GEO  to  MAG  is 


cos  0*cos(J)  cos  0*sm  (j)  -sin  0 

-sin  (J)  cos  <()  0 

= 

Uy 

sin  0«cos  <|)  sin  0»sin  cos  0 

rJ 

GEO 

N 

or 


^MAG  =  [(cos0cos(j))'V^  +  (cos0sin(t))-V^  +  (-sm0)'V^]  x  +  [(-sin4))-V^  +  (cos(j))-ig  y  + 
[(sin0cos(}))-V^  +  (sin0sin(t))'V'y  +  (cos0)-VJ  z 


The  parameters  of  the  centered  dipole  model  of  the  earth's  magnetic  field  specified  by  the  first  three  Gauss 
coefficients  gj®,  gj*,  h/.  The  formulas  for  determining  colatitude  and  longitude  of  the  intersection  of  the 
centered  dipole  axis  with  the  northern  hemisphere  are 


0 


and 


4>  =  tan'^ 


1*1 


The  values  of  colatitude  and  longitude  are  found  by  substituting  the  International  Geomagnetic  Reference 
Field  (IGRF)  or  the  Definitive  Geomagnetic  Reference  Field  (DGRF)  values  of  the  Gauss  coefficients  for 
the  desired  year.  These  values  can  be  found  in  Table  2  of  reference  12  and  Table  1  of  reference  13. 
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